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Peculiarities of length scales in a two-orbital superconductor 
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Abstract We study the spatial behaviour of coherency 
and magnetic field in a two-orbital superconductor. The 
superconducting phase transition is caused here by the 
on-site intra-orbital attractions (negative-?/ Hubbard 
model) and inter-orbital pair-transfer interaction. We 
find the critical (diverging at T c ) and non-critical (re- 
maining finite) coherence lengths and magnetic field 
penetration depth for various values of hopping inte- 
grals and the strengths of intra-orbital attractions. Nu- 
merical results have been obtained for a two-dimensional 
square lattice. 

Keywords Two-orbital superconductor ■ Negative-U 
Hubbard model ■ Coherence lengths • Magnetic field 
penetration depth 



1 Introduction 

The multi-component structure of electron system plays 
the crucial role in the formation of superconducting 
state in a number of novel compounds, see [TJ[2J[2J[3] . For 
the present time the multi-band nature of superconduc- 
tivity has been identified, among the others, in MgB2 
[5], cuprates [5], iron-arsenic compounds [7., VaSi [5], 
NbSe 2 0, and Sr 2 Ru0 4 PEHTT]. The presence of inter- 
acting superconductivity order parameters is decisive 
for these systems. 
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Recent studies [HHI3JII11IIS1II11III1IIHMI1] have re- 
vealed unexpected properties of two-band inhomoge- 
neous superconductivity compared to one-band case. 
Interband pair-transfer interaction modifies qualitatively 
the coherence lengths in two-band scenario p~3TfT4Tri6l 
ITT] . In this respect there is strong necessity of reinter- 
pretation of relevant experimental data. Moreover, the 
character of spatial coherency affects several theoreti- 
cal conception, e.g. type-1.5 superconductivity [13) sug- 
gested to explain novel patterns in vortex structure in 
magnesium diboride [20) and strontium ruthenate [21j . 

The characteristic lengths of a two-orbital super- 
conductor described by negative- U Hubbard model [23] 
were previously analyzed in dependence on band fill- 
ing for various inter-orbital interactions and orbital en- 
ergies p~8l!22J . As a continuation of that research we 
present below the results of the model calculations for 
the superconducting state and related length scales with 
varying strengths of intra-orbital interactions and hop- 
ping integrals. 



2 Basic equations 

The bulk superconducting state in a two-orbital system 
is described by the Hartree-Fock-Gorkov self-consistent 
equations 
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Here a = 1,2 is the orbital index; N is the number 
of lattice sites; A a are the superconductivity gaps; n a 
are the occupation numbers of orbitals per lattice site; 
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jjaa < q are ^Yie on-site intra-orbital attraction energies 
and U 12 = U 21 is the inter-orbital interaction energy. 
Note that both intra- and interorbital channels support 

superconductivity. Futher, E a (k) = \J £ 2 (k) + |^\ a | 2 , 

where s a (k) = e a (k) + \\J aa n a — fi and e Q (k) is the 
electron band energy associated with the orbital a; \i 
is the chemical potential determined by the equation 
n\ + n% = n for the total number of electrons per site 
n. The equations ([1]) and © follow in the mean-field 
approximation (see |18j ) from the Hamiltonian of the 
two-orbital negative-C/ Hubbard model 
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where is the electron creation operator in the or- 
bital a = 1, 2 localized at the site i\ a is the spin index; 
tf" is the hopping integral and e° a is the orbital energy. 
Eq. |T]) defines the superconductivity phase transition 
temperature T c by the equation 

(1 + U 11 g 1 (T c )) (1 + U 22 g 2 (T c )) 



(U 12 y gi (T c )g 2 (T c ) = 0, 
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For the spatially inhomogeneous two-orbital super- 
conductivity the Ginzburg-Landau equations 1181122] read 
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where for isotropic electron spectrum 

Pa = Pal = Pa2 = Pa3 

and 
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The expression for the density of supercurrent j s in- 
cludes the vector potential A, the magnetic flux quan- 
tum and the volume of unit cell Vq- 

The solutions of inhomogeneous gap equations © 
are scaled by two, critical or soft (index s) and non- 
critical or rigid (index r), characteristic lengths [TS] 



C,r (T) = 
with 

G(T) = ( 
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(6) Eq- © yields for the magnetic field penetration depth 



A(T) = 
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where /^o is the magnetic permeability of free space. We 
also consider here real homogeneous gaps A a (T). 

In what follows the numerical calculations have been 
carried out for two-dimensional square lattice with hop- 
ping integrals limited to nearest neighbours tf" = t a 
and orbital energies e± = e 2 = 0. In this case the elec- 
tron band energies associated with s-orbitals are e a (k) 
= — 2t a [cos(afc T ) + cos(aky)] where a is the lattice con- 
stant and —n/a < k x ,y < ix/a. We have also chosen 
k B = I- 
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Fig. 1 The critical temperature T c vs total number of elec- 
trons per lattice site n for the chosen set of hopping inte- 
grals (ti , t 2 ) and electron-electron interactions (U aa ; a, a' = 
1,2); (a) ti = t, t 2 = 0.8t, U 11 = -1.5t, U 22 = -2.5t, 
(b) ti = 0.8t, t 2 = t, U 11 = -1.5i, £/ 22 = -2.5i, (c) 
h = 0.8t, t 2 = t, U 11 = U 22 = -2.0t, (d) t t = t 2 = t, 
jyii = _ _ 2 .0t, respectively, and |t/ 12 | = |[/ 21 | = 0.04t 
in each case (t is an arbitrary unit of hopping integral). 




Fig. 2 Superconducting gaps A\ and A 2 vs temperature for 
the same sets of parameters as in Fig. 1 with the number of 
electrons n=2.5. The gaps coincide for the case (d) as the 
subbands (a = 1 and 2) are equal. 



3 Results and discussions 

In Fig. 1 the dependence T c vs n is shown for vari- 
ous values of U aa and t a . The shape of the function 
T c (n) is determined mainly by the following factors: (i) 
the location of the chemical potential fi in the separate 
bands including also Van Hove singularities (intraband 
contributions to superconductivity) and (ii) the loca- 
tion of /i in the region of bands overlapping (interband 
pair-transfer contribution to superconductivity). Note 
that the dependence T c (n) is symmetric for t\ = t% 
and U 11 = U 22 (curve (d) in Fig. 1). By introducing 
any difference in orbitals the asymmetry of the function 
T c (n) appears. The shift of the T c maximum towards 
the smaller values of n (curves (a), (b), and (c) in Fig. 
1) is caused by effective Hartree corrections U aa n a /2 
and by the difference in particular band- widths (8t a ) 
as U 11 U 22 and t\ ^ ti respectively. The stronger ef- 



fect of the Hartree renormalization |i[/ n ni 



is supported by the discrepancy of band fillings na,2 
which becomes especially large once /i passes the Van 
Hove singularity 

The temperature dependencies of the gaps A\ and 
A 2 (see Fig. 2) illustrate the variation of the driving role 



of pairing channels in different bands, including the tails 
(curves (a) and (b)) caused by the interband proxim- 
ity effect. In the latter case the intra-orbital processes 
are substantially more dominating in the formation of 
superconductivity compared to inter-orbital contribu- 
tion. In particular, there exists a temperature region 
just below the phase transition point where the pairing 
in the band with larger gap induces weak superconduc- 
tivity in the second band via interband pair-transfer of 
arbitrarily small intensity. 
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Fig. 3 The soft and rigid coherence lengths and penetration 
depth vs temperature. The parameters as in Fig. 1 with the 
number of electrons n=2.5. 
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The coherence lengths and magnetic field penetra- 
tion depth are plotted in Fig. 3. Apart from the diver- 
gence of the critical length scale £ s at T c , we observe also 
a maximum for the non-critical length scale £ r slightly 
below T c . This second peak is a combined effect of band 
filling and interactions in particular bands. 

It is worth to mention that for the parameters used 
there exists always a small temperature domain in su- 
perconducting phase where £ s > A > £ r . The domain 
width is determined by the values of A and £ r which 
are strongly affected by the interband proximity effect 
(see also Eq. (fTBjO and discrepancy between intrinsic 
critical temperatures of the bands. These factors play 
most in favour in case (a) . Note also that the presence of 
proximity effect can be identified by the non-monotonic 
temperature dependence of the critical coherence length 
£ s below T c (c.f. [17U18] h The corresponding maximum 
of £ s (T) for T < T c remains out off the scope of panels 
(a) and (b) in Fig. 3. 

4 Conclusion 

We have examined the temperature behaviour of crit- 
ical and non-critical coherence lengths, and magnetic 
field penetration depth in a two-orbital superconductor 
in dependence on orbital parameters. It was demon- 
strated that the variation of the combinations of intra- 
orbital electron-electron interactions and hopping in- 
tegrals changes substantially the relative efficiency of 
the pairing channels responsible for the superconduct- 
ing ordering which is reflected in the peculiarities of the 
dependencies of characteristic length scales on temper- 
ature. 
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